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Abstract — Solar cells have generated a lot of interest as 
a potential source of clean renewable energy for the future. 
However a big bottleneck in wide scale deployment of these 
energy sources remain the low efficiency of these conversion 
devices. Recently the use of nanostructures and the strategy of 
quantum confinement have been as a general approach towards 
better charge carrier generation and capture. In this article 
we have presented calculations on the optical characteristics 
of nanowires made out of silicon. Our calculations show these 
nanowires form excellent optoelectronic materials and may yield 
efficient photovoltaic devices. 

Index Terms — silicon, nanowire, solar cell, exciton 



I. Introduction 

IN recent years a lot of innovation in the semiconductor 
industry has enabled solar cell designs with higher 
conversion efficiencies though the use of inexpensive 
materials and processes [1], [2]. Today III-V semiconductor 
based solar cells can achieve upto 40% efficiency. Most 
strategies in improving the efficiencies target improving 
optical absorption, and efficient collection of generated 
photocarriers. In recent times improving efficiencies through 
quantum confinement has generated a lot of interest, with the 
advent of multijunction cells (MJSCs) and quantum dot based 
devices [3]. However most of the high efficiency devices 
are still based on exotic III -V hetrostructures, which are 
expensive to process and manufacture. Hence a real incentive 
in realizing cells made out of silicon still exists due to the 
existing infrastructure already in place. 

A possible candidate for the next generation of solar cells are 
silicon nanowires. Silicon nanowires are ID materials and 
there has been interest to incorporate them into photovoltaic 
cell design due to the wide availability of the material and 
existing processing infrastructure. The advantage of using 
nanowires being that, like quantum dots we can tune the 
band gap to a specific part of the solar spectrum, but with 
the added advantage that carrier transport in the solar cell is 
not dependent on quantum tunneling as proposed in [4]. In 
this paper we have presented the optical properties of silicon 
nanowires of different dimensions evaluated from ab initio 
calculations. Our results show that the nanowires are direct 
bandgap materials and we can alter the energy gap anywhere 
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between 2-5 eV, that gives us a wavelength selectivity from 
200 to 600 nm wavelengths. Results on intrinsic absorption 
and excitonic binding energies are also presented later in this 
article. 

The rest of this paper is organised as follows, Section-II 
describes the computational methodology followed for the 
ab initio calculations, In Section -III we present our results 
on bandstructure and optical properties, like the dielectric 
function and absorption coefficient. Finally we evaluate how 
exciton binding energies in these nanowires vary with the 
length of the nanowire (corresponds to quantum confinement 
in the axial direction). 



II. COMPUTATIONAL METHODOLOGY 

In this section we describe the first principles calculations 
carried out for the silicon nanowires systems. The 
computations were carried out using a density functional 
DFT-LDA approach using a parallel binary of the PWScf 
distribution. We carried out computations on three nanowires. 
We consider the nanowires to be oriented along the [100] 
plane, each having a square cross section. The dangling 
bonds on the surface of the nanowire are passivated by 
hydrogen terminations. Before performing electronic structure 
calculations we subjected the unit cells to a geometry 
optimization procedure. The equilibrium lattice relaxed to a 
unit cell with the Si-Si bond spacing to be 2.324 A and 
Si-H spacing to be 1.5 A. The translational lattice constant 
for the nanowires was 5.43 A. We have considered three 
nanowire unit cells in this paper each with a diameter 4.34 A, 
8.16 A and 10.75 A. These have been designated as Sill, 
Si22 and Si33 systems in this paper. The relaxed unit cells 
considered for the calculations presented in this paper are 
shown in Fig. 1. The silicon atoms in the Sill unit cell forms 
the primitive cell for the Si22, Si33 nanowires as well. The 
SiXX system is formed by 'X' repetitions of this atomic 
arrangement in the x and y directions. Each of the atoms 
in the Sill cell occupy four different z-planes at z = 0, 
z = ±ao/4 and z = a.o/2. 

The electronic structure of each of the nanowire unit cells was 
computed self consistently using a plane wave pseudopotential 
method. We utilized Van Barth Carr soft potentials for both 
silicon and hydrogen, with a Perdew Zunger exchange and 
correlation [5]. The SCF calculations were performed with a 
grid of 20 uniformly spaced k- vectors within the ID Brilliouin 
Zone of the nanowires. A kinetic energy cutoff of 80 Ry was 
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TABLE I: Computed Oscillator strengths for direct K=0 tran- 
sitions 
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Fig. 1: silicon nanowire unit cell geometries seen along 
the [100] direction. The three nanowires are designated (a) 
Sill,(b) Si22 and (c) Si33 in this paper 



deemed sufficient for numerical convergence of total energy. 
For the computation of optical absorption coefficient a(ui), a 
Drude-Lorentz broadening of 0.13 eV has been applied. 

Ill, Electronic and Optical Characteristics 

A. Bandstructure 

The silicon nanowires unlike bulk silicon are direct band gap 
semiconductors, which would make them excellent choices for 
optical applications. Authors in [6], [7], [8] have also reported 
direct band gaps and monotonically decreading bandgapes 
from their theoretical calculations. 

All the nanowires exhibit exceptionally high effective 
masses of holes and electrons. For example the Sill wire 
has a almost zero dispersion state on the top of the valence 
band. Similarly, the lowest conduction bands are also relatively 
heavy compared to their values in the bulk state. This is quite 
advantageous as it allows for extension of carrier lifetime 
as binding energy of excitons are directly dependent on it 
(5). We found that the energy gap of the nanowires decrease 
with increasing width. The wide tunability in the band gap 
of the nanowires allows us to tune the selectivity of optical 
absorption of the nanowire. The idea behind this being to 
match optical absorption in solar cells to the spectrum of 
the incoming solar radiation. A second intersting feature of 
the band structure is the weak dispersion of the conduction 
and valence bands leads to a high effective mass both in the 
conduction and valence bands. There is some experimental 
evidence [9] who have found out that carrier mobilities in in- 
trinsic nanowires is extremely low. Table-I shows the effective 
masses of the lowest four conduction bands CI — > CA for 
all the three systems. Due to the relatively weak dispersion 
seen in these wires especially the Sill and Si22 systems, it is 
prudent to present a thermally averaged effective mass, i.e. a 
harmonic average over a KT/q energy range from the lowest 
conduction band energy. 

B. Optical Absorption 

All the nanowires examined in this paper show features in the 
absorption spectrum the correspond to excitonic processes. A 
simultaneous comparison of their absorption spectra is shown 
in Fig. 4. The lowest excitonic peaks for each of the nanowires 
occur at 5.25 eV (232 nm), 3.7 eV (335 nm), and 2.3 eV (539 
nm) in the increasing order of size. The corresponding optical 
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wavelength is listed in parentheses. We see that absorption is 
tunable from the visible region to the near UV portion of the 
solar spectrum. These peaks are also quite prominent in the 
Sil 1 wire as compared to the Si33 nanowire due to the a small 
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Fig. 3: Optical Properties of silicon nanowires (a)-(c) Imaginary part of dielectric function 62(10) for the three nanowires shown 
for both longitudinal and transverse polarizations, (d) Comparison of refractive indices 
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Fig. 2: ID Bandstructure of silicon nanowire system along the 
z-direction with < k z < -k /arj with ag = 5.43 A being the 
lattice constant. Band diagrams are shown for (from left to 
right) Sill, Si22 and Si33 systems 



density of states in the conduction band. The rippled nature of 
the absorption spectrum weakens as the size increases. In the 
Si33 wire, these oscillations are almost washed out and the 
spectrum looks similar to that of a direct bandgap material 
with a exponential absoprtion edge near the conduction band 
energy [10], [11]. To ascertain the exact transitions that show 
up in the absorption spectrum, we computed the dipole matrix 
element of for each of the eigenvalues at k z = and computed 
a quantity called the oscillator strength /o. The quantity /o is 
proportional to the square of the dipole matrix as 



/o — -^tt (-^2 - Ei) 



ZTi 






(inlRM 



(1) 



where the summation is over all possible polarizations. 
The magnitude of /o is a direct gauge of the strength of 
absorption at a particular energy. These numbers have been 
tabulated in Table-I. We observed that only a few excitations 
have significant values of /rj appear optically active in the 
absorption spectrum. These excitations are also marked in 
Table-I. Other optical properties of these nanowires are also 
presented in Fig. 3. The strong transition peaks are more 
clearly visible in the plots for 62(0;). 
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Fig. 4: Optical Absorption coefficient for the nanowire systems 
a(ui) in units of energy (eV). 



C. Excitonic Binding energy 

It is known that excitonic processes are helpful in extending 
lifetimes of the electron hole pairs created by optical absorp- 
tion. The radiative lifetimes of an exciton r r is proportional to 

3/2 

the binding energy r r ex E b ' . Hence it is useful to evaluate 
how change in the length of the nanowire changes the binding 
energy of the ID exciton. To get this dependence we followed 
the methods of [12], [13] solving a two particle Schrodinger 
equation for an electron-hole system. In this section we have 
presented binding energy results for the lowest lying excitonic 
state in each of the nanowires. The Hamiltonian for such a 
system can be written in the Hartree approximation 



H = 



Pz 



n 2 
Pzh 



T 



(2) 
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The first two terms are simply the kinetic energy of the elec- 
tron and hole respectively and the third term is the screened 
Coulomb potential. The exciton is assumed to be immersed 
in a surrounding medium with an effective dielectric constant 
e r . For each of the nanowires we used a value computed from 
the real dielectric constant at zero frequency i.e. ei(0). This 
Hamiltonian is diagonalized by a two particle wave function 
of the form 



TTZ e 7TZ h _' 

ip = Jy cos cos e 



(3) 



where A is a fitting constant and N is a normalizing constant. 
This test wavefunction was originally applied to ID quantum 
wells. This trial wavefunction is acceptable in SiNWs as it 
is reasonable to conceptually think of a nanowire as a ID 
quantum well with an unsually long width. Equation (3) makes 
a single assumption that the wavefunction becomes zero at 
the ends of the quantum well, which is applicable to the 
nanowire as well. To find the binding energy of the exciton 
we use a variational method of minimizing the total energy 
over different values of A. The ground state of the exciton is 
simply 
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TABLE II: Radiative Lifetime of lowest optically active exci- 
tons 
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The minimized value of A corresponds to the effective Bohr 
radius of the exciton. The binding energies have been plotted 
as a function of nanowire length in Fig. 5. We see that the 
binding energies monotonically increase as the nanowire is 
made shorter. This dependence is expected as the stronger con- 
finement causes an increase in the effective mass of the charge 
carriers. The bulk excitonic binding energy is dependent on the 
reduced mass [i of the exciton given by the equation 



Ebulk = 



fiq 



(5) 



32ir 2 h 2 e 2 ,e 2 . 

where /i = m e mh/(rn e + m/,). We also see that the exci- 
tonic absorption is the strongest in the Sill system compared 
to the shorter excitonic peaks seen in the Si22 and Si33 system 
(Fig. 4). Using the method presented above, we computed 
the exciton binding energies for the strongest transitions in 
the nanowire systems. There are three strong transitions of 
interest in the Sill system and one transition each in the Si22 
and Si33 system (see Table-I). We found that the excitons 
in the nanowires are predominantly Mott Wannier type. The 
typical Bohr radii for the ID free exciton A/ was found to 
range between 2.5-5 nm. The largest radius being exhibited 
by the Si33 wire due to a strong screening potential. However 
confining the exciton in a short nanowire, with lengths less 
than Xf causes a large jump in the binding energies of the 
excitons as the electron hole pairs are forced to be closer to 
each other. The binding energies have a strong effect on the 
lifetime of these carriers. The radiative lifetime r is related to 
the oscillator strength and the exciton binding energy by the 
relation [14], [15] 
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where E v = E 9 — Eb, n is the refractive index of the 
material and the other symbols have their usual meanings. We 
found out that these nanowires exhibit unusually large exciton 
lifetimes with numbers of the order of a microsecond. We 
have summarized our findings for a specimen length of 10 nm 
(see Table-II). In the Sill wire designated Sillc the radiative 
lifetime is about 4 ns (I = 10 nm and similar vales of Eb), 
far less than the other two transitions Silla and Sillb. This is 
simply because of the large oscillator strength of this excitation 
is also responsible for a quick decay. 

IV. Conclusion 

We have presented first principles calculations for silicon 
nanowires for potential applications in photovoltaic applica- 
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Fig. 5: ID Exciton binding energy in silicon nanowires as 
a function of nanowire length. Excitons correspond to direct 
transitions labeled by a f in Table-I 



tions. Results indicate that the nanowires are direct band gap 
semiconductors and show excellent wavelength selectivity in 
the optical region of the solar spectrum. We also identified 
the optically active interband transitions present in these 
nanowires and used it to calculate excitonic effects in the 
nanowire by solving a two particle Schrodinger equation. Our 
calculations indicate that excitons in these nanowires show 
large radiative lifetimes of the order of several microseconds 
for nanowires as small as 10 nm in length. This number is 
expected to grow significantly with further confinement as we 
move into the quantum dot regime. 
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